The Distribution of Dark Matter in the Universe on scales of 10^*^ M(t, to 10^^ M.p 
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The use of paraUel computers and increasingly sophisticated software has allowed us to perform 
1 J - a large suite of N-body simulations using from 10* to 10^ particles. We will report on our recent 

r> I convergence tests of the halo mass function, N-point correlation functions, power spectrum and 

S ■ pairwise velocity from very large high resolution treecode N-body simulations. Rather than 

basing results on just one or two large simulations, now one can investigate the role of different 

numerical and physical effects on the statistics used to characterize the mass distribution of 

the Universe. 

1 Introduction 

Over the last two decades cosmological A^-body simulations have played a crucial role in the 
study of the formation and evolution of cosmic structure. In this approach, the initial statistical 
properties of the fluctuations, at some early stage of the universe, are set using linear theory. 
A'^-body simulations are then used to evolve the structure into the deeply nonlinear regime. This 
dynamical state is then compared with the large-scale structure in the galaxy data-sets. 
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The impressive progress achieved in the observational front with the completion of very large 
surveys such as 2dF and SDSS poses a clear challenge to the numerical work in cosmology: the 
precision of the predictions provided by the current A^-body experiments have to be of the order 
of a few percent. 

2 Numerical Simulations 

The cosmological model adopted is a low-density, flat, ACDM Universe with the parameters 
^m = 0.3, r^A = 0.7 and h = 0.7. The power soectrum of the initial conditions was set up 
using the output transfer function of CMBFASTI^, assuming $7;, = 0.04 and a normalisation 
of (Tg = 0.9. The entire suite of simulations were performed using the Hashed Oct-Tree code 
(HOT)'^, a parallel tree code with periodic boundary conditions. The simulation parameters 
are listed in Table 1. 

Table 1: Suite of simulations 
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3 Analysis 

The main purpose of this analysis is to understand how changes in mass resolution, particle 
number and box size affect some well known estimators applied in Large Scale Structure studies. 
For each cosmological volume we estimate the two-point correlation function, the power spectra, 
the mass function and the mean pairwise velocity, see Figure ^ The color map is shown in Table 
1. 

The two-point correlation function is presented in Figure ^(a). On small scales, the ampli- 
tude of the two-point correlation function is suppressed by softening gravity and lack of small 
scale power in the initial conditions. This affects the estimates in a range which roughly extends 
up to five times the size of the smoothing kernel [see arrows in top left panel of the abovemen- 
tioned figure]. The flattening of the two-point correlation function seems to disappear once the 
mass resolution (softening) increases (decreases). On scales where ^ ~ 1 all the simulations show 
a rather good agreement among themselves. The box size only looks to matter on scales well 
within the linear regime (^ -C 1) where the smaller boxes present a slightly smaller amplitude. 
The solid (dashed) lines represent the two-point Smith et al*^ (linear) correlation function. 



To measure the power spectra of our simulations we applied the technique detailed in Jenk- 
ins et al '^ and Smith et al '^. In doing so discreteness and grid effects are taken in account 
and removed from the final estimation. As expected, the power spectra of the different A^-body 
simulations show the same features as the two-point correlation function. The analytical ex- 
pression for the non-linear power spectra due to Smith et al*^ and the linear power spectra are 
shown by the solid and dashed black lines, respectively. Futher details can be found in Warren 
& Teodoro'"'. 

We used the friends-of-friends algorithm"^ to identify dark matter halos. This halo finder 
depends on just one parameter, b, which defines the linking length as bn~^'^ where n is the 
mean particle density; we follow the prescription of Jenkins et a/^and set b = 0.2. Figure (c) 
demonstrates that the agreement between the analytical expressions of Jenkins et al '^ (solid 
line) and Sheth and Tormenl^ (dotted line) and our estimates of the mass function is excellent. 
The small departures seen in the mass range M > 10^^/i~^MQ are mainly caused by small 
number statistics. 

Finally, Figuren(d) shows the mean pairwise velocity, vi2{r), as function of pair separation. 
The Hubble line, given by VHubbie = —Hr, is also shown. The mean pairwise velocity vanishes at 
the smallest separations resolved in our simulations. In the subset of simulations which better 
probe small scales (small volume and larger mass resolution) Vi2{r) closely follows the Hubble 
line up to ~ 300 kpc /i~^. This indicates that structures on such scales are close to relaxation. 
The mean pairwise velocity reaches a minimum in scales comparable to the correlation length [see 
Figure 1(d)]. At larger scales the vi2{r) intersects the Hubble line and, at very large separations 
it decays to zero, in accordance with the principle of large-scale homogeneity and isotropy. 

Some simulations in this suite have been further analyze in Seljak &; Warren"^. 
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Figure 1: (a) Two point correlation function, (b) power spectrum, (c) mass function and (d) mean pairwise 
velocity. See color map in Table 1. The arrows in the top left panel indicate the gravitational softening. In the 
same panel, we also show the linear theory (dashed line) and the Smith et al (solid line) two-point correlation 
function fitting for the assumed cosmological model. The solid, dotted and dashed lines in panel (c) show the 
Sheth-Tormen, Jenkins etal and Press-Schechter ^^ mass functions, respectively. Finally, in panel (d) the dashed 

line represents the Hubble flow. 



